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ABSTRACT 

We show that the recently-discovered orbits of massive stars with closest approach of 
~ 10 3 Schwarzschild radii from Sgr A* , allow winds from these stars to provide the 
required mass deposition rate near the black hole horizon. The observed luminosity 
of Sgr A* does not require viscous transport of angular momentum, as long as the 
total wind mass loss rate reaches ~ 1O _6 M0 yr _1 from the close-in stars. The specific 
orbits of the nearest stars should cause modulation of the radio and infrared flux from 
Sgr A* on a timescale of years in a predictable fashion. 
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1 INTRODUCTION 

Early in the development of accretion flow theory, it had 
been realized that angular momentum can severely limit 
the accretion rate of gas onto a black hole relative to the 
spherical infall value (the so-called Bondi (1952) rate). In- 
deed, most of the hot gas revealed by recent X-ray obser- 
vations within ~ 0.1 pc of the Galactic center (Baganoff ct 
al. 2003), is apparently not flowing into Sgr A* , the mas- 
sive black hole there (Quataert 2003). The latest models of 
the radiatively inefficient accretion flow around Sgr A* es- 
timate the accretion rate of gas into the black hole horizon 
at ~ 1O _8 M0 yr _1 , several orders of magnitude below the 
Bondi rate (Yuan, Narayan, Quataert 2003). This low ac- 
cretion rate is required by the frequency dependence of the 
linear polarization data (Bower et al. 2003), in combination 
with the synchrotron radio flux and the recently detected 
infrared emission from Sgr A* (Genzel et al 2003b; Ghez et 
al. 2003c). The absolute minimum accretion rate required 
by the bolometric luminosity of Sgr A* (~ 10 36 erg s _1 ) 
is ~ 2 x 10~ 10 M© yr~ 1 (e/0.1)~ 1 , where e is the radiative 
efficiency of the accreting mass. 

The disparity between the large gas reservoir available 
at the outer radius of influence of Sgr A* and the small mass 
that actually flows into its horizon, has inspired a variety 
of accretion models in which the mass inflow rate dimin- 
ishes with decreasing radius, due to convection (Quataert & 
Gruzinov 2000a, b; Narayan et al. 2000), turbulence (Coker 
& Melia 1997; Pen et al. 2003), or outflows (Begelman & 
Blandford 1999; Yuan et al. 2002). In this Letter, we high- 
light the significance of an important ingredient that was 
omitted by these models, namely the powerful winds from 



known stars that are embedded within the dilute accretion 
flow towards Sgr A*. 

Sgr A* was recently found to be surrounded by a cluster 
of young, massive stars. Over the past year, three of these 
stars, named SO-2, SO-16, and SO-19 by Ghez et al. (2003c; 
see parallel work by Schodel et al. 2003), were inferred to 
have orbits with peri-centric distances of ~ 100AU or ~ 10 3 
Schwarzschild radii from the central black hole. The detec- 
tion of a Bi'7 and He I features in of the spectrum of SO-2 
indicates that it is probably an O8-B0 star, while the other 
bright stars are also likely to be O-B stars, consistent with 
the lack of a CO absorption feature in their spectrum (Ghez 
et al. 2003a; Eisenhauer et al. 2003). 

The winds of O-stars on the main sequence have typical 
mass loss rates of M w — 1O _7 -1O _6 M0 yr _1 and speeds in 
the range v w — 1-3 x 10 3 km s _1 (Puis et al. 1996; Repolust, 
Puis, & Herrero 2003). In §2.1, we show that SO-2, SO-16, 
and SO-19 can supply the required low-angular momentum 
gas that fuels the emission from Sgr A* without requiring vis- 
cous transport of angular momentum. To demonstrate this 
point, we adopt the simplest model in which the radial gas 
streams move on ballistic Keplerian orbits. Gas with higher 
angular momentum is assumed, for simplicity, to be driven 
out of the system, as inferred to be case at the Bondi ra- 
dius (Quataert 2003). In §2.3 we show that the characteristic 
ram-pressure of the radial streams is well above the thermal 
pressure of the diffuse hot gas detected by the Chandra X- 
ray Observatory (Baganoff et al. 2003), making the ballistic 
approach viable. In §2.2 we illustrate how the specific orbital 
parameters of the close-in stars can be used to predict the 
long-term modulation of the radio and infrared flux of Sgr 



© 0000 RAS 



2 Loeb 



A*. Finally, we discuss the implications of these results in 
§3. 



2 DIRECT FLOW OF GAS INTO THE BLACK 
HOLE HORIZON 

2.1 Method of Calculation 

We consider the simplest model in which the gas elements 
move on ballistic Keplerian orbits. This assumption is justi- 
fied if: (i) viscosity is sufficiently weak so that mainly gas on 
nearly radial orbits approaches the horizon of Sgr A* ; (ii) 
the inflow of low angular momentum gas is dominated by one 
stellar wind at a time so that multiple gas streams do not 
collide; and (iii) the wind ram pressure is much higher than 
the gas into which it propagates, so that the low-angular mo- 
mentum gas is not deflected by hydrodynamic forces on its 
radial trajectory towards the black hole horizon. We justify 
the latter two conditions in §2.3. 

The mass of gas directed into the compact emission re- 
gion around the black hole horizon is the sum of all gas 
elements sent out from a stellar wind with angular momen- 
tum per unit mass smaller than some critical value, J m ax 
(see Fig. 1). We write, 




vsma 



M BH = if V 8> V » 



Figure 1. Schematic geometry of a radial stream. The stream 
is confined to a narrow cone centered on the direction where the 
tangential component of the wind velocity, v w sin a, cancels the 
tangential velocity of the star in its orbit around Sgr A* , vg (see 
Eq. 2). The small opening angle of the radial stream, S, is dictated 
by the maximum angular momentum, J m ax, that allows the gas 
to join the compact emission region around Sgr A* (Eq. 3). 



(1) 



where M = 4 x 1O 6 M is the black hole mass (Ghez et 
al. 2003b; Schodel et al. 2003), and AGM/c is the max- 
imum angular momentum with which the gas would en- 
ter directly into the horizon of a Schwarzschild black hole 
(see equation f4.2.f2 in Shapiro & Teukolsky f983). For a 
nearly radial trajectory starting at the star, we find that 
f] = 0.5x/(x ~ f) 1/2 , where x = r min /(2GM/c 2 ) is the dis- 
tance of closest approach to Sgr A* (in Schwarzschild radii 
units) of the ballistic gas orbit*. In popular emission models 
(Melia et al. 200f; Melia, Liu, & Coker 200f; Liu & Melia 
2002a; Yuan et al. 2002; Yuan et al. 2003), most of the 
bolometric luminosity of Sgr A* is produced within ~ 30 
Schwarzschild radii, implying r\ <> 2.8 for the wind gas to 
join the emission region directly from its orbit. Next, let us 
calculate the fraction of the wind material that has an or- 
bital angular momentum \ J\ < J m ax for a star at a radius r 
in an orbit with an eccentricity e and a semi-major axis a. 

The low angular-momentum gas is confined to a narrow 
cone centered on an angle a relative to the direction of Sgr 
A* (see Fig. 1), such that 



v e + v w sin a = 0, 



(2) 



where v w is the wind speed, ve — ijr is the tangen- 
tial component of the orbital velocity of the star (using 
the polar coordinates (r, 6) in the orbital plane), and I = 



[GMa(l - e 2 



11/2 • 



is the orbital angular momentum of the 



* This result is obtained by setting E 1 and solving V(r) 
in equation 12.4.25 of Shapiro & Teukolsky (1983). 



star. Note that condition (2) may be satisfied only through 
a portion of the stellar orbit that is sufficiently far from Sgr 
A* so that v w is larger than vg oc r _1 ; elsewhere, we set the 
mass deposition rate into Sgr A* to zero (see the gaps of 
vanishing Mbh in Fig. 2). The small opening angle of the 
Sgr A* acceptance cone, 5 <S a, is defined by the condition, 



J ir 



ve + v w sin(a + 8) 

Using the small-angle approximation, sin(a + 8) 
8 cos a, we find 

c- i/max 



(3) 



sin a + 



v w r cos a 



(4) 



where cos a = [l — (£/v w rY 



,1/2 



We can now evaluate the fraction of the wind mass that 
is channeled through the acceptance cone of the emission 
region around the black hole horizon. Assuming that the 
wind outflow is spherically-symmetric in the star rest-frame, 
this fraction is 



Mr 



_ / Jmax \ 

47T _ \2v w r) 



(5) 



For the characteristic values of v u 



1000 km s" 1 and r 



lO- 3 AU, one gets M BH ~5x W- 1 (-n/3yM w . 

If the gas element leaves the star at a time t, then it 
will reach the black hole on a radial orbit at a time t + At. 
The time delay At can be easily calculated from the initial 
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Figure 2. Mass deposition rate into the emission region sur- 
rounding the horizon of the black hole at the galactic center, 
Mbh> as a function of calender year. Different line types cor- 
respond to the individual contributions of the three bright stars 
nearest to Sgr A* , namely SO-2 (solid), SO-16 (dashed), and SO- 
19 (dotted). The orbital parameters of these stars were adopted 
from Table 3 of Ghez et al. (2003). The sensitivity of the re- 
sults to the assumed wind velocity, v w , is illustrated by the 
difference between the top (v w = 1000 km s -1 ) and bottom 
(v w = 3000 km s" 1 ) panels. The gaps in Mbh occur during 
phases with v$ > v w along the stellar orbit, and are centered 
on pericenter passages. 



radial distance of the star r, and the initial radial velocity 
of the gas, 



dr 

U r = —± V u 

at 



-i 1/2 



(6) 



where there are either zero, one or two allowed cones under 
the restriction that u r < 0. The timescale for the gas to 
be drained into the black hole horizon from the compact 
emission region around it (Hawley, Balbus, & Stone 2001), 
is expected to be short compared to the timescales of interest 
here and so we neglect it. 



2.2 Numerical Results 

Figure 2 shows the numerical solution of equations (5) and 
(6) for Mbh as a function of t + At, given the orbital pa- 
rameters of SO-2 (solid line), SO-16 (dashed), and SO-19 
(dotted), as summarized in Table 3 of Ghez et al. (2003). 
The evolution of the mass deposition rate into the black 
hole horizon is presented for the values v w = 1 x 10 3 km s" 1 
(upper panel) and 3x 10 3 km s _1 (lower panel), that bracket 
the range of wind speeds deduced from observed spectra of 
O stars (Puis et al. 1996; Repolust et al. 2003). 

For the value of M w ~ 10" 6 M yr _1 found in bright 
Galactic O stars on the main sequence (Puis et al. 1996; 
Repolust et al. 2003), we get a total feeding rate of ~ 



10 (??/3) Mq yr , comparable to the value required by 
the polarization data in accretion models of Sgr A* (Yuan 
et al. 2003; Melia, Liu & Coker 2000; Bromley, Melia, & 
Liu 2001). Most of the wind material spreads over the 
highly eccentric orbit of its parent star, and so we find that 
M w <, 10" 6 M @ yr _1 for SO-2, SO-16, or SO-19 yields elec- 
tron densities below the canonical inflow model of Yuan et 
al. (2003) with the resultant Faraday rotation measure sat- 
isfying current polarization constraints. 

It is natural to expect the luminosity of Sgr A* to scale 
with Mg H , where /3 ranges between 1 (for disk or jet emission 
with a constant radiative efficiency) and 2 (for radiatively 
inefficient flows). The emission spectrum can be generically 
explained by synchrotron and inverse-Compton processes 
(Melia & Falcke 2001; Yuan et al. 2002; Yuan et al. 2003) in a 
hot, magnetized gas within tens of Schwarzschild radii from 
Sgr A* , independent of the exact accretion geometry. The 
radial streams that reach this region would likely dissipate 
a fraction of their kinetic energy through shocks, and settle 
into a rotationally supported configuration around Sgr A* 
(see Hawley, Balbus, & Stone 2001, for the subsequent evo- 
lution of such as torus) . The interaction of the freely- falling 
material with the rotating torus may resemble the situa- 
tion in wind-fed X-ray binaries (Beloborodev & Illarionov 
2001). The magnetic fields carried by the infalling streams 
could be amplified by plasma instabilities and produce the 
synchrotron luminosity of Sgr A*. 



2.3 Are Ballistic Orbits Plausible? 

The properties of the accretion flow under consideration here 
resemble the conditions in wind-fed X-ray binaries (Illari- 
onov & Beloborodev 2001, and references therein). The ac- 
creting matter in such binaries originates in a wind from an 
OB star with similar mass loss rates to those assumed here, 
and is thought to freely-fall towards the boundary of the ac- 
cretor (be it the surface of a neutron star or the horizon of 
a black hole) where the X-rays are produced. The ballistic 
approach adopted here had been used successfully to derive 
the accretion rate and the corresponding X-ray luminosity 
of accreting neutron stars (Illarionov & Beloborodev 2001). 

Figure 2 illustrates that at most times there is only one 
gas stream dominating the mass deposition rate into Sgr 
A*. Collisions of multiple streams with comparable momen- 
tum flux are rare. The dominant gas streams follow ballistic 
orbits as long as their ram pressure, p w — p w v^, exceeds 
the thermal pressure of the dilute gas through which they 
propagate, p t h- We find, 



M w v w 3 • 

P-w = -; — 7T = 3 x 10 M w - 6 v w ,2 



where 



10 3 AU 



kcV 



(7) 



M w ,- 6 = (M w /10- b M e yr" 1 ), 

v m> 2 = («u>/2000 km s" 1 ), and d is the distance from the 
star (ignoring gravitational focusing by Sgr A* ). The dif- 
fuse X-ray emitting gas, spread over a scale of Jj 10 4 AU 
from Sgr A* , has a characteristic temperature of ~ 1.3 keV 
and a density ~ 30 cm -3 (Baganoff et al. 2003), implying 
a background thermal pressure, p t h ~ 40 keV cm" 3 , which 
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is below the ram pressure of the radially infalling streams 
at r <§C 10 4 AU. As each stream falls towards the black hole, 
it gets accelerated and compressed and so its ram pressure 
increases. Since its infall speed scales as oc r~ 1//2 and its 
surface area must decrease as r~ 2 , we expect that the ram 
pressure would scale as p w oc r~ 3 and that the condition 
Pw > Pth will be maintained throughout the trajectory of 
the stream. Our simplifying assumption (which we intend 
to critically examine with a direct numerical simulation in 
the future) is that most of the outer gas detected by Chan- 
dra or any other gas with significant angular momentum is 
not able to flow inwards similarly to the radial streams in 
the absence of efficient angular momentum transport. Note 
that a clumpy wind would enhance the survivability of the 
radial streams and potentially lead to short-term variability 
of Sgr A*. 

The fate of wind gas with angular momentum above 
Jmax remains an interesting open question. The linear po- 
larization data (Bower et al. 2003) implies that most of this 
gas does not reach the vicinity of the black hole. Since most 
of the energy dissipation occurs at ~ 10 Schwarzschild radii 
and the gravitational binding energy per particle decreases 
considerably at larger radii, we conjecture that heat flux 
from the inner regions unbinds the outer gas and expels it 
in an outflow (similarly to the conditions at the Bondi ra- 
dius; see Quataert 2003). 



3 DISCUSSION 

We have found that radial gas streams from the winds of 
massive stars with recently discovered orbits that approach 
within ~ 10 3 Schwarzschild radii of the black hole at the 
Galactic center, can supply the required mass flow near 
the central black hole horizon and power the radio and in- 
frared emission of Sgr A*. For O star winds with a total 
M w ~ 1O _6 M yr" 1 , we get that SO-2, SO-16, and SO-19 
can supply as much as ~ 1O _8 (?7/3) 2 M0 yr _1 of gas in orbits 
directed straight into the emission region around the black 
hole horizon. The total mass inflow rate would increase if an- 
gular momentum of other gas elements is transported away 
by magnetic or turbulent stresses. The surprising existence 
of massive stars in orbits with closest approach distances of 
~ 100 AU poses a major puzzle to theories of star formation 
(Genzel et al. 2003a), and was not incorporated in previous 
modeling of the gas flow around Sgr A* (e.g. Coker & Melia 
1997; Yuan et al. 2003). 

The parameters of the winds could in principle be de- 
termined through detailed spectroscopic observations of the 
stars in the Sgr A* cluster. Once the wind properties of the 
individual stars are constrained, our accretion model will be 
easily testable since it correlates the long-term (>month) 
variability of the radio or infrared emission from Sgr A* with 
the orbital phases of the stars around it (see Fig. 1). The 
simplified calculation presented in this Letter can be im- 
proved using a full hydrodynamic simulation and compared 
against observational data on the long-term variability of 
Sgr A* (Zhao et al. 2003) for different choices of v w and M w 
per individual star (for an alternative source of this variabil- 
ity, see Liu & Melia 2002b). A proper numerical treatment 



is challenging since only a small fraction of the gas mass 
follows the radial orbits of interest here, and since the in- 
flow spans several orders of magnitude in radius - requiring a 
large dynamic range in spatial and temporal resolution. Pre- 
vious accretion simulations (such as those by Coker, Melia, 
& Falcke 1997, Proga & Begelman 2003 or Rockefeller et 
al. 2003) did not meet the demanding requirements of this 
accretion problem. 

Sgr A* may be surrounded by fainter stars that have 
not been detected so far, even though they dominate the 
steady mass flow into the black hole. However, winds from 
massive stars are driven by radiation pressure, and so the 
stars dominating the fluctuations in the mass flow into 
the black hole horizon should have the highest luminosi- 
ties and be included in the current census of bright in- 
frared sources near Sgr A*. Massive stars that are much 
further away also possess powerful winds (Krabbe et al. 
1991; Najarro et al. 1997; Coker & Melia 1997; Quataert 
2003), but their individual contributions to the radial mass 
flux near the horizon of Sgr A* is reduced compared to 
the close-in stars^, since Mbh oc r~ 2 in equation (5). For 
example, IRS 16C with its M w w 8 x 1O~ 5 M yr _1 and 
v w ~ 650 km s" 1 at r ^ 2 x 10 4 AU, and IRS 13E1 with 
its M w « 8 x 10~ 4 M Q yr _1 and v w ~ 1000 km s _1 at 
r ^ 4 x 10 4 AU (Quataert, Narayan, & Reid 1999), pro- 
duce at most a radial mass flux comparable to that of SO- 
2 if the latter has M w ~ 10~ 6 M Q yr -1 . Nevertheless, the 
steady mass flux into Sgr A* may still be dominated by un- 
detected stars or by gas that originates at large distances 
and loses its angular momentum through viscous transport 
(Yuan et al. 2003). In this case, the winds from the close- 
in stars would only modulate a fraction of Mbh, but their 
fractional contribution could still be identified through the 
predictable short-term variability that their orbits introduce 
to the luminosity of Sgr A*. 
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